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TITLE OF THE INVENTION 
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CLAIM OF PRIORITY 

The present invention claims priority from 
Japanese application JP 2003-114508 filed on April 18, 
2003, the content of which is hereby incorporated by 
reference on to this application. 

BACKGROUND OF THE INVENTION 

The present invention relates to a biomagnetic 
measurement technique detecting a very weak magnetic 
field such as a magnetocardiogram and a 
magnetoencephalogram using an SQUID (superconducting 
quantum interference device) magnetometer as a 
superconducting device . 

An action potential occurs by excitation of a 
ventricular muscle to produce a very weak electric 
current in the ventricular muscle in appearance. Along 
with the electric current, a very weak magnetic field* 
is produced outside a living body. The very weak 
magnetic field from a heart (hereinafter, abbreviated 
as a ''magnetocardiogram'') is measured using a highly 
sensitive magnetic sensor called a superconducting 
quantum interference device (SQUID) . The 

magnetocardiogram is less affected by the conductivity 
of a living body, which is hard to be subject to 
waveform distortion. From the measurement result of the 



magnetocardiogram, an attempt to image an electric 
current activity on the surface of a ventricular muscle 
has been made. 

In the attempt to image an electric current 
activity, there is proposed a method for calculating Ix 
= dBz/dy and ly = -dBz/dx from a magnetic field (Bz) of 
a vertical (z) element to the surface of a living body 
to reconstruct a distribution of imaginary current 
vectors I = (Ix, ly) (for example, see Japanese Patent 
Application Laid-Open No. 10-248821) . 

As an ischemic area identification method, 
there is proposed a method for calculating the total of 
absolute values of the imaginary current vectors I in a 
predetermined period (for example, see Japanese Patent 
Application Laid-Open Nos 10-305019 and 11-151221).. An 
ischemic area is estimated by mapping of the total 
value in the predetermined period. 

As an invasive method for directly measuring a 
ventricular muscle action potential, there is catheter 
mapping. The catheter mapping is a method for directly 
contacting a catheter having at its edge a plurality of 
potential measurement electrodes with the inner surface 
of a heart to measure a ventricular muscle action 
potential, 

SUMMARY OF THE INVENTION 

In the above prior art, the amount of electric 
current can be calculated by a distribution of current 



vectors and the total of absolute values of the current 
vectors in a predetermined time. No waveforms 
corresponding to a ventricular muscle action potential 
can be measured. The catheter mapping in the prior art 
is a method for inserting a catheter into a body under 
X-ray radiation, which is an examination putting a 
large load on patients. 

An object of the present invention is to 
provide a biomagnetic measurement technique which can 
obtain a potential waveform corresponding to a 
ventricular muscle cell action potential in a non- 
invasive manner. 

To achieve the above object, in a biomagnetic 
measurement apparatus of the present invention, when a 
plane in parallel with a plane contacted with the 
surface of a living body is xy plane and an axis 
vertical to the xy plane is z, a magnetic field 
produced from a heart, is detected by a plurality of 
magnetometers including superconducting quantum 
interference devices . The superconducting quantum 
interference devices are driven by ah operating circ^uit 
for magnetometer. Output data (signal data) of the 
operating circuit for magnetometer are collected by a 
data collection device. An operation processor executes 
operation processing of the collected output data 
(signal data) . The result of the operation processing 
is displayed on a display device. 

Ix = dBz/dy and ly = -dBz/dx are calculated 



from a detected magnetic field (Bz (x, y) ) of a 
vertical (z) element to the surface of a living body to 
calculate imaginary current vector I = (Ix, ly) . An 
absolute value of the current vector is calculated by 
using Ixy (t) = V"{ (Ix (t))^ + (ly {t))^}. 

When a potential waveform at start time to of 
depolarization of the heart (a period in which QRS wave 
appears) is V (to) = 0 or Ixy (to) = 0 and end time of 
repolarization is tn, a potential waveform at time ti (i 
= 0, 1, . . . , m) in a period from the to to end time tm 
of depolarization of the heart is calculated by using V 
(ti) = V (t i-i) + Ixy (ti) . The value of potential 
waveform Vo (tm) at the end time tm of depolarization is 

A potential waveform at the start time t ni+i of 
repolarization of the heart (a period in which ST-T 
wave appears) is calculated by using V (t m+i) = - Ixy 
(t m+i) . A potential waveform at time ti (i = m+2 , 
m+3, ..w n) in a period from the time t ni+2 in a period 
of repolarization to the end time tn of repolarization 
is calculated by using V (ti) = V (t i-i) - Ixy (ti) . 

The method for obtaining the imaginary current 
vector I is not limited to calculation of Ix = dBz/dy 
and ly = -dBz/dx. A current vector calculated by a 
solving method of inverse problem using a lead field 
i^^trix and a current vector using aminimum norm method 
may be used. 

More specifically, the operation processor 



executes the following first to third operation 
processing. 

In the first operation processing, when t = ti 
(i 0, 1, . . . , m) is a period corresponding to 
depolarization of the heart 'of the living body' and t = 
ti (i = m+1, m+2 , n) is a period corresponding to 

repolarization of the heart of the living body-v from 
the element Bz in the z direction at time t at the 
measurement point (x, y) , a current vector (Ix (t) , ly 
(t) ) at time t and an absolute value of the current 
vector (Ixy (t) (Ix (t) ) ^ + (ly (t))^} ) is 

calculated at time ti (i = o, 1, n) . 

In the second operation processing, when the 
lower limit of addition E is i = 0 and the upper limit 
of addition Sisi.=0, 1, ...,m, a potential waveform 
at time t = ti (i = 0, 1, m) in a period 

corresponding to depolarization of the heart of the 
living body is calculated by using V (ti) = S Ixy (ti) . * 

In the third operation processing, when the 
value of the potential waveform V (tm) at the end of a 
period corresponding to depolarization of the heart of 
the living body is Vm, the lower limit of addition E is 
i = m+1, and the upper limit of addition E is i = m+1, 
m+2, . . .', n, a potential waveform at time t = ti (i = 
m+1, m+2, . . . , n) in a period corresponding to 
repolarization of the heart of the living body is 
calculated by using V (ti) = Vn, - E Ixy (ti) . 

The' potential waveform V (x, y, t) at time t = 



ti (i = 0, 1, n) at the measurement point (x, y) , 

which is obtained in the second and third operation 
processing, is displayed on the display device. The 
operation processor calculates the potential waveforms 
5- V (x, y, t) at a plurality of measurement points (x, y) . 
The equipotential diagram connecting the equipotential 
points of the potential waveforms V (x, y, t) ' 
calculated at the plurality of measurement points (x, 
y) is displayed by contour map. The waveform of the 

10 measured magnetic field and the potential waveform V (x, 
y, t) are displayed on the display device. 

The biomagnetic measurement apparatus and the 
data processing method of the present invention can 
obtain an action potential waveform corresponding to 

15 each area of a heart. Information on an abnormal 

ventricular muscle action potential such as long QT 
syndrome can be obtained in a non-invasive manner. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram showing the configuration 
example of a biomagnetic measurement apparatus of an 
embodiment of the present invention; 

FIG. 2 is a diagram showing an example of an 
array of detection coils arranged in the cryostat of 
25 FIG. 1; ' 

FIG. 3 is a diagram showing an example of 
magnetocardiograms measured by the biomagnetic 
measurement apparatus of FIG. 1; 
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FIG. 4 is a diagram showing a. current vector 
calculation method of the present invention; 

FIG. 5 is a diagram showing the relation 
between a cell-membrane ionic current, a ventricular 
muscle action potential and an EGG according to the 
embodiment of the present invention; . 

FIG. 6 is an equivalent circuit diagram of 
assistance in explaining a mechanism about electrical 
propagation of adjacent ventricular muscle cells 
according to the embodiment of the present invention; 

FIG. 7 is a diagram schematically showing an 
example of a method for calculating a ventricular 
muscle action potential waveform according to the 
embodiment of the present invention; 

FIG. 8 is a diagram of assistance in explaining 
an example of calculation processes calculating a 
ventricular muscle action potential waveform according 
to the embodiment of the present invention; 

FIG. 9 is a diagram showing an example of 
ventricular muscle action potential waveforms 
calculated from the magnetocardiogram waveforms of FIG. 
3 ; and 

FIG. 10 is a diagram showing comparison of a 
ventricular muscle potential measurement result (the 
top figure) by a catheter examination measured in the 
right ventricular wall of a patient having Type I long 
QT syndrome, an overlapped figure (the middle figure) 
of the magnetocardiogram waveforms, and the calculated 



action potential waveform (the bottom figure) according 
to the embodiment of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
A biomagnetic measurement apparatus of an 
embodiment of the present invention in which when a 
plane in parallel with a plane contacted with the 
surface of a chest is xy plane and an axis vertical to 
the (x, y) plane is z, current vectors from magnetic 
fields produced from a heart and absolute values of the 
current vectors are calculated at a plurality of 
measurement points (x, y) , wherein potential waveforms 
at the plurality of measurement points (x, y) are 
calculated based on current vectors in a period from 
start time of depolarization of the heart to end time 
of depolarization of the heart and absolute values of 
the current vectors in a period from the end of 
depolarization of the heart to end time of 
repolarization of the heart, thereby detecting 
information on early af terdepolarization of the heart. 

In a period from start time of depolarization 
of the heart to end time of depolarization of the heart, 
absolute values of the current vectors are added to 
calculate a potential waveform. In a period from the 
end of depolarization of the heart to end time of 
repolarization of the heart, the absolute values of the 
current vectors are subtracted from the value of the 
potential waveform at the end of depolarization of the 



8 



heart to calculate potential waveforms at the plurality 
of measurement points (x, y) . 

More specifically, when start time of 
depolarization of the heart is to, end time of 
depolarization of the heart is tm, end time of 
repolarization of the heart is tn, an absolute value of 
a current vector at time ti is Ixy (ti) , and a potential 
waveform at the time to is V (to) = 0, a potential 
waveform at time ti in a period from the time to to time 
tm is calculated by using V .(ti) = V (t i-i) + Ixy (ti) . 
The value of potential waveform V (t„,) at the time t^ is 
Vra, and a potential waveform at the time t n»+i is 
calculated by using V (t m+i) = - Ixy (t m+i) . A 
potential waveform at time ti in a period from time t ^+2 
to time tn is calculated by using V (ti) = V (t i-i) - 
Ixy (ti) . The potential waveforms V (ti) are calculated 
at the respective measurement points (x, y) . 
Information on early af terdepolarization can be 
detected from the potential waveforms V (x, y, t) 
calculated at the plurality of measurement points (x, 
y) . The thus-calculated potential waveform V (x, y, ti) 
corresponding to each area of the heart is displayed on 
the display device. 

A data processing method in the biomagnetlc 
measurement apparatus of the present invention having a 
plurality of detection coils detecting, when a plane in 
parallel with a plane contacted with a living body is 
xy plane and an axis vertical to the (x, y) plane is z, 



element Bz in z direction of a magnetic field produced 
from the living body, and a data collection device 
collecting detected signal data of the element Bz in 
the z direction, which is executed after collecting the 
signal data, wherein in a period from start time of 
depo-larization of the heart to end time of 
depolarization of the heart, absolute values of current 
vectors are added to calculate a potential waveform, 
and in a period from the end of depolarization of the 
heart to end time of repolarization of the heart, the 
absolute values of the current vectors are subtracted 
from the value of the potential waveform at the end of 
depolarization of the heart to calculate potential 
waveforms at a plurality of measurement points (x, y) . 

More specifically, when start time of 
depolarization of the heart is to. end time of 
depolarization of the heart is t^, end time of 
repolarization of the heart is tn, an absolute value of 
a current vector at time ti is Ixy (ti) (i = 1, 2, 
n) , and a potential waveform at the time to is V (to) = 
0, a potential waveform at time ti in a period from the 
time to to the time is calculated by using V (ti) = V 
(t i-i) + Ixy (ti) . The value of potential waveform V 
(tm) at the time tm is Vn,, and a potential waveform at 
the time t is calculated by using V (t m+i) = - Ixy 
(t m+i) . A potential waveform at time ti in a period from 
the time t n.+2 to the time tn is calculated by using V 
(ti) = V (t i-i) - Ixy (ti) . The potential waveforms V 
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(ti) are calculated at the respective measurement points 
(X, y) . 

An embodiment of the present invention will be 
described below in detail based on the drawings. 

FIG. 1 is a diagram showing the configuration 
example of a biomagnetic measurement apparatus 
according to an embodiment of the present invention. In 
a magnetically shielded room 1, there are arranged a 
bed 7 on which a living body is placed, a cryostat 2 
storing a refrigerant (liquid helium or liquid 
nitrogen) for holding SQUID sensors in a 
superconducting state, and a gantry 3 for fixing the 
position of the cryostat 2. The SQUID sensors are 
operated as magnetometers by a driving circuit 4 
arranged outside the magnetic shield room 1. The 
outputs of the magnetometers pass through an amp filter 
unit 5 to be converted to digital data by an A/D 
converter circuit incorporated in a computer 6, and are 
then stored in the computer 6. 

FIG. 2 is a diagram showing an example of an 
array of detection coils 8 arranged in the cryostat 2 
of FIG. 1. The detection coils 8 are integrated with 
the SQUID sensors. In the example shown in FIG. 2, the 
detection coils 8-1 to 8-64 integrated with the SQUID 
sensors are arranged in an 8 X 8 matrix. The arrangement 
in a matrix facilitates calculation of a current vector 
of the present invention (which will be described later 
in detail) . As the current vector calculation method. 



there are a minimum norm method and a method for 
calculating a lead field inverse matrix. The detectiSdn 
coils are not necessarily arranged as shown" in FIG. 2. 

FIG. 3 is a diagram showing an example of a 
magnetic field waveform (hereinafter, called a 
^'magnetocardiogram waveform") measured in the heart of 
a patient having Type I long QT syndrome by the 
biomagnetic measurement apparatus of FIG. 1, The top 
figure of FIG. 3 shows a grid map 9 showing magnetic 
field waveforms corresponding to the respective 
positions (channels) of the detection coils 8-1 to 8-64 
shown in FIG. 2. The bottom figure of FIG. 3 shows an 
enlarged figure 10 of the grid map in four positions in 
which characteristic magnetic field waveforms appear. 
From the enlarged figure 10, in the lower parts of the 
detection coils in four positions (channels) , it is 
found that two-layer (notch type) characteristic 
waveforms appear. In FIG. 3, the double-headed arrow 
line on the horizontal axis indicates a time width of 1 
sec, and the double-headed arrow line on the vertical 
axis indicates 40pT. 

FIG. 4 is a diagram showing a process creating 
a current vector from the z element (Bz (x, y) ) of a 
magnetic field used' in the embodiment of the present 
invention. The top figure of FIG. 4 shows an overlapped 
waveform 11 of magnetic field waveforms in which the 
magnetocardiogram waveforms (the grid map 9) of 64 
channels shown in FIG. 3 are overlapped with each other 



on one trace to be displayed. 

Using a distribution 13 of Bz (x, y) at an 
observation time 12 of a distribution 15 of current 
vectors In (x, y) , calculation 14 of the current vector 
In is performed by Ix = dBz/dy and ly = ~dBz/dx. This 
can obtain the distribution 15 of current vectors In (x 
y) showing directions and magnitudes by straight lines 
with arrows. The distribution 15 of current vectors In 
shows actual data about a patient having Type I long QT 
syndrome. It is found that abnormal current (the black 
part in the figure) is produced in the lower part of 
the measured plane of the magnetocardiogram. In FIG. 4, 
the double-headed arrow line on the horizontal axis 
indicates a time width of 1 sec, and the double-headed 
arrow line on the vertical axis indicates 50pT. 

As the current vector calculation method, the 
minimum norm method and the method for calculating a 
lead field inverse matrix may be used. 

FIG. 5- is a diagram showing the relation 
between a cell-membrane *ionic current, a ventricular 
muscle action potential and an EGG according to the 
embodiment of the present invention. The distributions 
of ion concentrations are different inside and outside 
the ventricular muscle cell. When the permeability of a 
cell-membrane ion is selectively advanced, the ion is 
flowed in or out of the cell according to the 
electrochemical potential gradient. The ion flowed in 
and out of the cell membrane can be considered as an 



electric current, which is called a cell-membrane ionic 
current. The entrance of the cell-membrane ionic 
current is called an ion channel. FIG. 5 shows a 
representative cell-membrane ionic current. 

In depolarization period A, an inward current 
(an electric current flowing from outside the cell to 
inside the cell) 16 having a large amount of electric 
current is flowed inside the cell in a short time. As 
the representative inward current 16, there is Na + 
electric current (INa) . After that, the ventricular 
muscle is brought to plateau phase (refractory period) 
B during which period there is no significant cell- 
membrane ionic current movement. Finally, it is brought 
to repolarization period C so that an outward current 
(an electric current flowing out from inside the cell 
to outside the cell) 17 is slowly flowed out from 
inside the cell. As the representative cell-membrane 
ionic current of the outward current 17, there is + 
electric current (IK) . It is found that a large number 
of ions are involved in the inward current 16 and the- 
outward current 17. FIG. 5 is a schematic diagram 
showing the cell-membrane ionic current very simply. 

A ventricular muscle action potential 18 is 
formed by the cell-membrane ionic currents (16, 17) 
flowed out from outside the cell or to inside the cell. 
The ventricular muscle action potential 18 is a 
waveform performing synthesis in the present invention. 
The bottom figure. of FIG. 5 shows an EGG 



waveform 19 finally measured as the potential of the 
surface of the living body by the total electric 
activity of the ventricular muscle action potential 18. 
It is considered that the ECG waveform is formed by the 
differential potential between the action potential of 
the endocardium and the action potential of the 
epicardium of a ventricular muscle cell or the 
differential potential (electrical propagation) between 
the action potentials of adjacent ventricular muscle 
cells. As a result, in the ECG waveform 19, QRS 
waveform appears in the depolarization period A and T 
wave appears in the repolarization period C / 

FIG. 6 is an equivalent circuit diagram of 
assistance in explaining a mechanism about electrical 
propagation of adjacent ventricular muscle cells . 
according to the embodiment of the present invention. 

FIG. 6 shows three-dimensional ventricular 
muscle tissue as a two-dimensional propagation model 
with three ventricular muscle cells 20-1, 20-2 and 20-3 
Only electrical propagation at depolarization will be 
described here. When an electric signal transmitted 
from the stimulation transmission system of the heart 
in the left direction is transmitted to Vq and Vo in the 
ventricular muscle cell 20-1 exceeds the threshold 
value (about -60 to -70mV) , the switch of an ion 
channel 21-1 is opened to start flowing cell-membrane 
ionic current imi as the inward current via a parallel 
circuit having resistance r^i and capacitance 0^1, The 



flowing of the cell-membrane ionic current imi produces 
an action potential specific to the ventricular muscle 
cell 20-1 in. the Vq. 

Subsequently, electric current iii2 is flowed 
via junction resistance rua to a gap junction part 22-a 
of the ventricular muscle cells to produce an action 
potential in VI in the ventricular muscle cell 20-2. As 
in the ventricular muscle cell 20-1, when the VI 
exceeds the threshold value (about -60 to -70mV) , the 
switch of an ion channel 21-2 is opened to start 
flowing cell-membrane ionic current 1^2 as the inward 
current via a parallel circuit having resistance and 
capacitance Cni2 . The cell-membrane ionic current im2 
produces an action potential specific to the 
ventricular muscle cell 20-2 in the Vi . 

Subsequently, electric current ii23 is flowed 
via junction resistance ri23 to a gap junction part 22-b 
of the ventricular muscle cells to produce an action 
potential in V2 in the ventricular muscle cell 20-3. The 
excitation is propagated to the adjacent ventricular 
muscle cells. It is supposed here that signals trapped 
in the magnetocardiogram are the electric currents iii2 
and ii23 considered to be the propagation signals. 

Under the supposition, the electric currents in, 
the ventricular muscle at a certain time are considered 
to be iii2 = (Vi - Vo)/rii2 and ii23 = (V2 - Vi)/ri23, which 
reflect the potential difference between the cells. To 
calculate an action potential of each ventricular 



muscle cell, an action potential at time at which cell 
excitation does not occur must be zero to calculate the 
total amount of electric current corresponding to 
ventricular muscle electrical excitation sequentially 
produced. It is considered that a potential 
corresponding, to an action potential can be calculated 
by the total amount of electric current. In the above 
description, due to calculation at depolarization, the 
total amount of electric current is considered as the 
action potential calculation method. At repolarization 
for the outward current (in the direction opposite the 
inward current at depolarization) , there is performed a 
method for subtracting from the total of electric 
currents (which will be described in detail in FIGS. 7 
and 8) . 

FIG. 7 is a diagram schematically showing a 
method for obtaining a ventricular muscle action 
potential waveform according to the embodiment of the 
present invention. FIG. 7 shows a schematic diagram of 
the action potential calculation described in FIG. 6. 
FIG. 6 assumes that a current vector is calculated in a 
certain magnetocardiogram measurement position 
(channel) . In the depolarization period (QRS waveform) 
A of a magnetocardiogram waveform 23, as shown in an 
enlarged figure 25 of the obtained action potential 
waveform, absolute values of current vectors are added. 
The absolute values of the current vectors are added in 
the + direction to obtain a total value. 



In time from the start of the refractory period 
B to the end of the repolarization period C, the 
absolute values of the current vectors are subtracted 
from the largest value 2 6 at the end of the 
depolarization period (QRS waveform) A. The absolute 
values of the current vectors are subtracted (or added 
in the - direction). As a result, an action potential 
waveform 24 is obtained. 

FIG. 8 is a diagram of assistance in explaining 
calculation processes obtaining a ventricular muscle 
action potential waveform according to the embodiment 
of the present invention. 

In the following description, t = ti (i = 0, 
1, . . . , m) is a period corresponding to depolarization 
of the heart of the living body and t = ti . (i = m+1, 
m+2, n) is a period including the refractory 

period B to the end of repolarization of the heart. 

From a magnetic field measured at time t = ti (i 
= 0, 1, n) , a current vector (Ix (t) , ly (t) ) and 

an absolute value of the current vector (Ixy (t) = ^^ 
{ (Ix (t)).^ + (ly (t))^} are calculated (a process 33 
calculating an absolute value of a current vector) 

In a period of depolarization of the heart, 
absolute values of current vectors are added. When the 
lower limit of addition S is i = 0 and the upper limit 
of addition 2 is i = 0 , 1 , . . . , m, a potential 
waveform at time t = ti (i = 0, 1, . . . , m) in a period 
corresponding to depolarization of the heart is 
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calculated by using V (ti) = E Ixy (ti) . Here, the 
junction resistance of the gap junction part is 
considered to be 1. 

More specifically, V (to).= Ixy (to); V (ti) = 
Ixy (to) + Ixy (ti) ; V (tz) = Ixy (to) + Ixy (ti) + Ixy 
(t2) ; V (t„,) =.Ixy (to)) + Ixy (ti) + ... + Ixy (t^^) 

= Vm are calculated (a first process 34 calculating 
potential wavef orm -V) . is the largest value 2 6 of the 
potential waveform at the end of the period 
corresponding to depolarization of the heart. 

In a period from the start of the refractory 
period B to the end of the repolarization period C, the 
absolute values of the current vectors are subtracted 
from the largest value 26 at the end of the 
depolarization period (QRS waveform) A. When the lower 
limit of addition E is i = m+1 and the upper limit of 
addition E is i = m-hl , *m+2 , n, the absolute 

values of the current vectors are subtracted to 
calculate a potential waveform at time t = ti (i = m+1, 
m+2, n) by using V (ti) = Vm - E Ixy (ti) . 

More specifically, V (t m+i) = V xy (t ni+i) ; V 

(t ni+2) ) = - Ixy (t m+l) " Ixy (t m+2) ; . . . ; V (tn) = 

Vra - Ixy (t mti) - Ixy (t zn+2) " ..." Ixy (tn) are 
calculated (a second process 35 calculating potential 
waveform V) . 

In the result of the second process 35, a 
display process 36 performs waveform display of the 
potential waveform V (the bottom figure of FIG. 9) , 



grid map display (the top figure of FIG. 9) , and 
contour map display of an equipotential diagram 
connecting the equipotential points of the potential 
waveform V. Data for diagnosis are provided as display 
easily understood by an examiner. Data about the 
contour map display are omitted. 

FIG. 9 is a diagram showing ventricular muscle 
action potential waveforms obtained from the 
magnetocardiogram waveforms of FIG. 3. FIG. 9 displays 
action potential waveforms 27 about 64 channels 
calculated using actually measured data of the patient 
having Type I long QT syndrome shown in FIG. 3. In 
particular, the action potential waveform of the area 
in which the action potential waveform 27 shows a 
characteristic pattern is shown as an enlarged figure 
28. The action potential waveform showing a 
characteristic pattern appears in an area almost 
corresponding to a right ventricle. 

According to the embodiment of the present 
invention, FIG. 10 shows, on one trace at the same time, 
a ventricular muscle potential measurement result 29 of 
one channel by a catheter examination measured in the 
right ventricular wall of the same patient (having Type 
I long QT syndrome) with the measured magnetic field 
waveforms (magnetocardiogram waveforms) of FIG. 3, a 
magnetocardiogram waveform overlapped figure 30, and an 
obtained action potential waveform 31 (the potential 
waveform in the right ventricle shown in the enlarged 



figure of FIG. 9) . Time on- the horizontal axis is 
indicated by standardized QTc having RR interval of 1. 
One or more magnetocardiogram waveforms and one or more 
ventricular action potential waveforms are displayed on 
the same slide at the same time to easily understand 
the relation between the magnetocardiogram waveform and 
the ventricular action potential waveform. 

When observing the ventricular muscle potential 
measurement result 29, a notch can be found in the 
calculated action potential waveform 31 at time at 
which early af terdepolarization (EAD) 32 occurs. The 
ventricular muscle potential measurement result 29 
coincides well with the calculated action potential 
waveform 31. They are found to . have been measured in 
the same right ventricle. The effectiveness of the 
ventricular muscle action potential waveform 
calculation method according to the embodiment of the 
present invention can be understood. The present 
invention can obtain information on an abnormal 
ventricular muscle action potential in a non-invasive 
manner without performing a catheter examination. 

The present invention can calculate an action 
potential waveform corresponding to each area of a 
heart and obtain information on an abnormal ventricular 
muscle action potential in a non-invasive manner. 



